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ABSTRACT 

We study the dynamical evolution of a stellar disk orbiting a massive black hole. We explore 
the role of two-body relaxation, mass segregation, stellar evolution and binary heating in af¬ 
fecting the disk evolution, and consider the impact of the nuclear cluster structure and the 
stellar-disk mass-function. We use analytic arguments and numerical calculations, and apply 
them to study the evolution of a stellar disk (similar to that observed in the Galactic center; 
GC), both on the short (few Myr) and longer (100 Myr) evolutionary timescales. We find the 
dominant processes affecting the disk evolution are two-body relaxation and mass segregation 
where as binary heating have only a little contribution. Massive stars play a dominant role in 
kinematically heating low mass stars, and driving them to high eccentricities/inclinations. 
Multi-mass models with realistic mass-functions for the disk stars show the disk structure to 
be mass stratified, with the most massive stars residing in thinner structures. Stellar evolution 
plays an important role in decreasing the number of massive stars with time, thereby leading 
to slower relaxation, where the remnant compact objects of these stars are excited to higher 
eccentricities/inclinations. At these later evolutionary stages dynamical heating by the nuclear 
cluster plays a progressively more important role. We conclude that the high eccentricities and 
high inclination observed for the majority on the young O-stars in the Galactic Center sug¬ 
gest that the disk stars had been formed with initially high eccentricities, or that collective 
or secular processes (not explored here) dominate the disk evolution. The latter processes are 
less likely to produce mass stratification in the disk; detailed study of the mass-dependent 
kinematic properties of the disk stars could therefore provide a handle on the processes that 
dominate its evolution. Finally, we find that the disk structure is expected to keep its co¬ 
herency, and be observed as a relatively thin disk even after 100 Myrs; two-body relaxation is 
too inefficient for the disk to assimilate into the nuclear cluster on such timescales. It therefore 
suggests earlier disks now containing only older, lower mass stars might still be observed in 
the Galactic center, unless destroyed/smeared by other non-two-body relaxation processes. 

Key words: Galaxy; centre - galaxies: nuclei - Galaxy; structure - stars; kinematics and 
dynamics - Galaxy: nucleus 


1 INTRODUCTION 

Stellar disks are known to exist around massive black holes 
(MBHs) in galactic nuclei, and are typically composed of relatively 
young stellar populations. Observations of the central parsec of the 
Milky-Way Galaxy s how more than 100 young WR, O and B stars 
jPaumard et alj200^ . likely formed ~ 4 — 7 Myr ago. These stars 
are found between ~ 0.04 — 0.5 pc from the MBH, with ~ 20% of 
which reside in a disk-like structure, and are observed to have mean 
high eccentricity of (e) ~ 0.3 — 0.4. Stellar disks might therefore 
be a frequent phenomena in galactic nuclei. Such disks are thought 
to form following the infall of gaseous material, and its formation 
of a gaseous disk. Such a disk c an later fragment to form a stellar 
disk jHobbs & Navakshinll2009l . and references therein). Here we 


explore the dynamical evolution of such stellar disks, and focus on 
disks similar to that observed in the Galactic center. 

[Alexander et alj ( l2007h have been the first to describe the be¬ 
havior of a stellar disk around a MBH. They used an analytic ap¬ 
proach, and then verified and calibrated them through the use of N- 
body simulations. The analytic model they devised includes the ef¬ 
fects of two-body relaxation, and included simple multi-mass mod¬ 
els (two or three populations). Their N-body simulations also in¬ 
cluded a simplified stellar evolution. 

Later studie s bvICuadra et al.l ( l2008h . |Perets et aP jlOOSl) and 
iLockmann et alj ( I2OO9I) performed additional N-body simulations 
and c onsidered additiona l components and/or more realistic as¬ 
pects. ICuadra^T^ ([200^ included stellar binaries in the disk, the 
role of a putative intermediate mass black holes, and initially ec- 
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centric disks. iKocsis & Tremaiii^ ( 1201 ih discussed the evolution 
of the st ellar disk, mainly focusing on the effects of resonant re¬ 
laxation dRauch & Tremainall996li . not considered h ere. The in¬ 
fluenc e of the stellar cusp was stud ied numerically bv IPerets et al.l 
(l2008h and lLdckmann et ^ l l2009h . who performed N-body simu¬ 
lation for different models (including models with cusp and models 
with two stellar disks); they concluded that the cusp could provide 
an additional heating of the stellar disk, leading to faster eccentric¬ 
ity/inclination growth. 

Here we follow -up on the initial analytic work by 
[Alexander et al.l (l2007h . and extend it to include the role of a re¬ 
alistic detailed mass-function, the effects of stellar evolution, the 
impact of binary-heating and additional heating by a stellar cusp. 
This approach allows us to better understand and identify the role 
played by each of these processes and components, as well as to 
study the long-term evolution of realistic disks, which is more dif¬ 
ficult (and computationally expensive) to study using full N-body 
simulations. 

We begin by describing the dynamical processes that con¬ 
tribute to the stellar disk evolution, including stellar scattering, 
mass segregation and binary-heating (section 2). We then briefly 
review the structural components of the nuclear stellar cluster in 
which the disk is embedded, and consider various models for the 
mass function of the stellar disk, and the possible core-like or 
cuspy structure of the nuclear stellar cluster (NSC). The physical 
processes we consider as well as the NSC structural components 
are used to devise a multitude of possible models for the stellar 
disk evolution under diverse conditions, including considerations 
of stellar evolution (section |4j. We present the results of the short 
and long term evolution of the disks considered in these models in 
section [3 and then discuss the results and their implications (sec¬ 
tion!^. We also briefly mention other possible physical processes 
not included in our current modeling. Finally we conclude and sum¬ 
marize our results in section|7] 


2 THE DYNAMICS AND PHYSICAL PROCESSES IN A 
STELLAR DISK 

In the following we briefly review several dynamic al processes that 
contrib ute to the disk evolution. We first follow [Alexander et al.l 
([2007[) in describing the two-body relaxation and mass-segregation 
(dynamical friction) of multiple stellar populations (extended to an 
arbitrary number of stellar populations), and then consider the an¬ 
alytic description of the effects of binary heating on stellar disks, 
not discussed before in this context. 


2.1 Relaxation and single-star scattering 

Let us follow [Alexander et al.[ l [2007h and consider a super massive 
black hole (MBH) with mass M, and N single stars with mass M 
orbiting the MBH in a initially thin stellar disk, with an external 
radius Ro + AR/2 and an internal radius Ro — AR/2. 

Wh en M, ';$> NM the relaxa tion time of such a system is 
given bv [Binnev & Tremaine! ([2008[) . 


-C^2A^pln(A) 

With p the stellar density and C ~ 1 (depending only 
on the system geometry; e.g. for a spherical case C ~ 0.34 
( [Binnev & Tremain3[2008[) . The disk density is given by 



Figure 1. Evolution of the velocity dispersion of stars in a stellar disk 
compo sed of single mass star's orbiting a MBH; similar to [Alexander et al] 
<2007li . 


NM 

2tiRoHAR' 


( 2 ) 


where H = a/Qis the scale height of the disk. We then obtain 


^ ^ Ro2ARo'^ 

trela:, Q2 ^ In (A) ’ 


(3) 


Using this Eq. and defining torb = 2 ti/Q. we can describe the 
evolution of the velocity dispersion of disk stars: 


da G^iVM^ In (A) 

dt CitorbRoARaS’ ^ ’ 

with Cl = 2C. 

In a disk of stars orbiting an MBH, the root mean square 
(RMS) ecc entricity of the stars ca n be related to their velocity dis¬ 
persion by ([Alexander et al.ll2007[ . and references therein) 

( c ) — firms — ^^2 , 

VK 

where vk = {GM,/Ro)^^^ is the Keplerian orbital speed. The 
RMS inclination of stars is expected to be of the order of half the 
RMS eccentricity value. For comparison with observations of the 
eccentricities of stars in the Galactic center we therefore present 
our results of the RMS eccentricity (and inclination in some cases) 
rather than the velocity dispersion. 

Fig.dlshows the evolution of the RMS eccen tricity of stars in 
the ste llar disk, using a similar setup as used bv [Alexander et all 
( [2007[) 


2.2 Multi-mass components and mass segregation 

Stellar scattering had also been considered for multi-mass stellar 
popul ations in the dispersion dominated regime (|Alexanderrt_^ 
[200^ . In this regime, for two mass populations. [Goldreich et ah 
( 20041) suggest that the treatment of the problem only depends on 
the amplitude of the velocity dispersion ( 0 - 2 ) of the lighter objects, 
with mass M 2 , compared with the Hill velocity of the more massive 
components , with mass Mi. The Hill velocity for the more massive 
stars is defined as 
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Vh,i = ^Rh,i, (5) 

where Rh,i = Ro{Mi/M 2 f^^■ The dispersion dominated regime 
occurs when a 2 > vh,i ■ In this regime the speed of the single stars 
is similar to the velocity dispersion and scattering encounters are 
well approximated hy two body dynamics. In the shear dominated 
regime, not discussed here, it is necessary to take into account the 
tidal gravity of the central hlack hole. 

When the velocity dispersion of the low-mass stars satisfies 
Vsac,i > 0-2 > vh,i, with Uesc.i the escape velocity from massive 
stars, the exchange of momentum between the heavier and lighter 
bodies occurs through collisionless gravitational deflectio ns, where 
gravitational focusing is important 1 Goldreich et al.ll2003l . 

Considering only two stellar populations composed of Ni and 
N 2 stars with stellar masses Mi a nd M 2 , and velocity di spersions 
(Ti and (T 2 , respectively, we obtain ( [Alexander et al.ll2007h : 

dai NiMi A^ 2 iW'iM 2 In (A 12 ) (Ji / 

~dr “ Altorh^rl Miff, W2\ MJ 

( 6 ) 

d(j 2 _ N 2 M 2 -1 A'iMiM 2 In (A 12 ) cr 2 A 

dt Altorb<^i ^ A2torb ^^2 V E 2 J ' 

(V) 

where cti 2 = (cri -|- (T 2)/2 and Ai = {CiRoAR)/G^ for i = 1 
and 2, and Ai A 2 are the appropriate Coulomb logarithms for the 
two stellar populations,which depend on the respective scale height 
of the disk. A 12 is related to the scale heig ht at which interactions 
between the two stellar species occur (see [Alexander et ^[20071 
for further details) The constant C 2 is related to the constant Ci, 
and also depends on the geometry of the system, i.e for a three 
dime nsional system C 1 IC 2 — 3.5 (for a detailed analysis please 
see at [Stewart & Wetherill[[l98^ : here we adopt the same value. 

We can extend such a model to an arbitrary number of dif¬ 
ferent stellar species, N. Each species adds an additional coupled 
equation and contributes a coupling term to each of the other equa¬ 
tions. For the general case the time evolution of the velocity dis¬ 
persions for each of the stellar species is given by a set of coupled 
equations 


daj NjMj ^ NkMjMkln{Ajk) Oj / Ek\ 

dt ~ AlGrbO^^ ^ A2torb V Ej) 

( 8 ) 

for i = 1, 2,. .., and k — 1,2,... ,N (with k 7 ^ j). 

From analyzing these set of equations we note that for two 
species, k and j for which Ek > Ej, the population of fc—species 
particles dynamically heat the population of the y-species ones 
(leading to increased velocity dispersion) and vice versa (a less en¬ 
ergetic population dynamically cools a more energetic one, dump¬ 
ing its velocity dispersion). 


stars interact with those cluster stars which go through the disk, i.e 
if the volume of the disk is 1/ = AttRqH AR , the number of clus¬ 
ters stars considered is the nuclear cluster stellar density times this 
volume. The velocity d ispersion of this isotropic c luster population 
is (T cluster Vk / teinnev & Tremain3[2008h with Vk the ke- 
plerian velocity around the MBH. We now use the same type of 
two-population interaction term in Eq.[ 8 ]to account for the effects 
of the cluster stars on the dis stellar population. Note that the clus¬ 
ter stellar population is much greater than the disk population, and 
can be assumed as a constant “heat” bath. 


2.4 Binary Heating 

Any stellar population is likely to contain some fraction of binary 
stars. An encounter between single stars and binaries could lead 
to an energy exchange between the binary orbital energy and the 
kinetic energy of the stars, i.e., energy can be exchanged between 
the inner degree of freedom of the binary and the outer external ki¬ 
netic energy of the stars; thereby binaries can assist in dynamically 
heating the stellar disk. 

Let us consider a binary system composed of stars with masses 
mi and m 2 , center of mass velocity I4m , and orbital velocity V 
with respect to the center of mass of the system (we assume a cir¬ 
cular binary for simplicity). The total e nergy of the binary can be 
expressed as teinnev & Tremain3[2008l) : 


E=-yy-G 


mim2 , mi-I-m2, ^2 

“T y nr 


where n — is the reduced mass. The sum of the first two 

^ mi+m 2 

left element define the binding energy, Eb, in terms of the semi¬ 
major axis 


Eb = -G 


mim 2 

2 a 


(9) 


Binary-single encou nters had been exp lored extensi vely, with 
the pioneering studies bv [Heggi3 l [l975ll an d [Hills[ l [l975h . Binaries 
can be divided into two dynamical categories; “Hard Binaries” and 
“Soft Binaries”: when the kinetic energy of the intruding mass, Ek, 
is such that EklrUbV'^ <C 1 , (where mt, = mi -|- m 2 is the binary 
total mass)a binary is considered to be hard (where the lower the 
ratio the harder the binary); it is considered to be a soft binary when 
Ek / 

[Heggie[ | [l975l) and [Hillsl ( [l975h found that following a binary- 
single encounter soft binaries become softer and hard binaries be¬ 
come harder (on average) , a behavior typical termed as “Heggie’s 
law”. 

The long term typical evolution of a binary randomly encoun¬ 
tering stars in some dense environment can then be modeled using 
a simple approach. Considering a binary in a field of incoming stars 
with typical mass m and particle density n , the averaged evolution 
of the binding energy of the binary due to many encounters is given 
by 


2.3 Interaction between the stellar disk and the nuclear 
cluster 

The two-body relaxation of disk stars due to interaction with the 
background stars of the nuclear cluster can be modeled in a similar 
method to the two-body relaxation by different stellar populations 
discussed above. In this case, the cluster stellar population is con¬ 
sidered as a distinct population with different properties. The disk 


=n{aEV)Eb, ( 10 ) 

where V the typical relative velocity of the intruders, and as the 
cross section of the single-binary collisions. We may now con¬ 
sider the contribution of a stellar population i with stellar density 
m which is composed of stars with typical mass rm, to the hi- 
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nary evolution; when the velocity distribution of in coming sta rs is 
Maxwellian and taking Eq.[T0]can be expressed as ( lHillslll992h : 

g-g (67r)2 G^mlpi f rm + f 1 \ 

V rn, ) + 

( 11 ) 

where mb = mi + m 2 , pi = rumi. iHillsI lll992h made use 
of N-Body simulations to study the dependence of on the in¬ 
truder mass (mi) and the binary mass (mb) in the mass-ratio range 
2mijmb € [0.01,10000] and found that the the empirical func¬ 
tional dependence can be formulated by 

* _ a + cx + ex^ -f gx^ 

^ 1 -I- fea; -I- dx^ -I- fx^ 

with X = \n{2mi/mb), a = 1.26, b = —0.15, c = 0.054, 
d = 0.02, e = 0.1, / = 0.0045 and g = 0.015. 

For a single star the average kinetic energy is (E) = Sma^ /2 
with (Ji the velocity dispersion. Taking ct; ~ 1 / (V^) (encounter 
velocity ^dispersion velocity) and applying energy conservation 
for a system of binary stars with density ribwe get 



rii 




Eqs. [TT] and [T^ give 


( 12 ) 


with 



rib 

Zmmiaf 


(13) 


Di = 



/ mj -h mb 

V rub 



We can now use Eq.[T^as the binary-heating term to be added 
to the dynamical evolution equation for the stellar disk populations 
(Eq.ll; allowing us to model the binary-heating contribution. 

We are now in a position to asses whether binary heating can 
play an important role in heating the stellar disk.. For the simplest 
case of a disk with only one population of same-mass stars and 
a binary fraction rib, the velocity dispersion evolution equation is 
obtained by summing the scattering term and the binary heating to 
get (taking mb = mi = M) 


da _ G^NM'^ In (A) rib 

m ~ CitorbRoARa^ imMa^’ 

or after some math 


^ _ G^NM^ nn(A) 
dt torbRo^Ra^ Cl ^ N 
where Nb is the number of binary stars, and 


(15) 


7 

We can now compare the relative contribution by single-star 
scattering and binary heating through the comparison of the two 
terms in Eq. [15] If we consider some typical initial conditions in 
the Galactic center, taking In (A) ~ 7.5 (as the disk heats-up the 


scale height grows but with only a small corresponding change in 
ln(A) due to the logarithmic dependence) , Ci ~ 3 and a binary 
fraction of 2/3, Nb = (2/3)W, we find that the relative contribu¬ 
tion of binary heating is only one tenth of the contribution due to 
single star scattering. This can be seen explicitly in the evolution 
of the velocity dispersion shown in Fig. The left panel shows 
how the exclusion of the binary heating contribution (for a disk in 
which Maingie. = Mbinary = 20Mq) makes only a very small 
difference (see also the effect of accounting only for binary heat¬ 
ing). We can therefore conclude that although the binaries binding 
energy could be larger that the kinetic energy of the disk, the rate 
of extraction of this energy and its exchange into the disk kinetic 
energy is negligible, at least for the relatively thin disks considered 
here. 


2.5 Full model 

We can combine the stellar disk contributions described above to 
obtain the full set of equations describing our model, including 
multi-mass populations, binaries and nuclear cluster stellar popu¬ 
lation. This is summarized below; 


dt 


NjMf 


r^itorba^ 


In (Aj 


NkMjMkln{Ajk) aj 


A2torb 

X ^ pv '^binary— i 


jk 


1 - 


Ek 

Ei 


(16) 


-E 


^binary— i-^^jT^binary — i Ih — binary — i) 


j — binary — i 


1 - 


-^ 2 ^orb 
Ebinary — i 


NcuspMjMcusp In (A j — cusp') 


A2to 


^ Ecusp 

Ei 


where 


(i) Aj—cusp — aj—cuBpARf2G]A[Q 

(ii) aj—cusp — (acusp -f rrj)j2, 

(hi) acusp = I i/GMt (R) /R 

(iv) Mt (R) = Mcusp (R) + Msmbh 

(v) Mcusp (R) = 47rr^pcusp (r) dr , where we con¬ 

sider the effective region to affect the stars to be between Rmin = 
0.5R andRmax = 1.57? and 

(vi) Pcusp (r) oc 

/ ^.7 —orbitg/ 3 

(vii) Aj —binary —i ~ ‘^GMy^nary — i ’ ^^Cre Vi —orbital ~ 

-and Oi is semi-major axis of the binary. 


(viii) Ajk 


a /^R 
2 GMfc 


3 STRUCTURAL COMPONENTS OF THE NUCLEAR 
STELLAR CLUSTER 

In the following we describe the observed properties of the clus¬ 
ter and the embedded stellar disk, and the models we use in our 
analysis. 
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Figure 2. Evolution of RMS eccentricity for model X (in Section |4| and 
Table[3]for more information about the models). . In the top panel is plotted 
the evolution of RMS eccentricity through time, the dashed line is the model 
only considering auto-relaxation and the solid one is when binary heating 
is included. In the bottom panel are ploted the the evolution of the RMS 
eccentricity when only auto-relaxation (dashed line) or binary heating is 
considered (solid line). 

3.1 Observed disk structure and mass function 

Beginning with the analysis of iLevin & Beloborodov! l l2003h . ob¬ 
servations have shown the exis tence of a young stel lar disk near 
the MBH in the Galactic center. Il^umard et aP j2006h showed that 
many of the observed young OAVR/B-stars reside in a disk-like 
structure between ~ 0.05 — 0.5 pc from the MBH. 

Early measurements suggested th e disk-stars revolve around 
the MBH in somewhat eccentric orbits iB^tko et al] ( l2009l) found 
an RMS eccentricity of (e) — 0.37 ± 0.07 for the disk stars, but 
the most recent stu dies suggest a smaller value of (e) = 0.27±0.07 
jYelda et alj2014 ). The disk-stars appe ar to be young, forming only 

4 — 7 Myrs ago < Paumard et al.ll2OO0 . The radial surface density 
profile of the disk can be described by a po wer-law distribution , 

5 (R) ^ R~^ , with 0 in the range [1.5 , 2.3] l lPaumard et al .120061 : 
iLuetal .1120091: iBartko et al.ll20od 1201 Oh . 

Initially it was thought that ~ 40 — 50% of the young stars 
belong in the disk structure; however, more recent studies suggest 
that only ~ 20% of these reside in the disk, while most of the star s 
have more spherical distribution, outside the disk jLu et aljl2009h . 
Here we consider only the evolution of disk stars. 

In most galactic environments the initial mass function (IMF) 


describes the di stribution of stellar masses at birth. Analysis 
of observations dSalnetej Il955l . IChabriej 1200 ll . iKroupal 1 20021) 
shows an almost universal power law behavior, with dN (m) oc 
m-“( with a typical exponent of a ~ 2.3, for stars above 

one Solar mass an d up to ~ 120 Mq. The analysis done by 
iBartko et al.l ( l2010h based on spectroscopic survey made with SIN- 
FONI conclude that the present day mass function and the IMF for 
the disk of young stars in the Galactic center obey a power law, but 
the slope is flatter than the “regular” IMF. In particular, they con¬ 
clude that a = 0.45 ± 0.3. In such a “top heavy” IMF, massive 
stars are much more frequent compared with the Salpeter IMF. 

3.2 Disk models 

In our models we consider an initially thin cool disk {H -C Ai? 
and low initial eccentricities for the star, i.e. we assume an initially 
low velocity dispersion for the disk stars) rotating at angular speed 
= VkIR o (with Vk the Keplerian velocity). 

We only consider the evolution of the disk stars and do not 
discuss the origin and evolution of the more spherical component 
of the young stars in the Galactic center. 

We studied three different types of disks: a single-mass disk, 
an Initial Mass Function (IMF) disk and disk with binaries. 

3.2.1 Mass function 

We consider various mass functions for the stars in the stellar disk 
and the cusp. These different models are summarized in the follow¬ 
ing. 

Single-mass model: In the simplest, unrealistic model studied 
here, we consider a stellar disk compo sed of single-mass stars , sim¬ 
ilar to the basic models considered bv lAlexander et aU l l2007t) . as a 
test-case for comparison. 

Top-heavy initial mass-function model: In this model we assume 
the initial mass function of the disk stars is a top-heavy mass func- 
tion with masses i n the range 0.6 — 120 M©, as suggested by 
IBartko et al.l 1 I 2 OIOL see discussion above). 

Salpeter mass function: In this case we assume the initial mass 
function of the stars in the stellar disk is a Salpeter mass function 
with masses in the range 0.6 — 120 Mq . 

3.2.2 Binarity 

Our knowledge about binaries in th e GC is currently li mited, but 
Observation by lOtt et al.lll99^ and iMartins et alj [20061 show the 
existence of at least some young binaries in the Galactic center. 
In particular, they found that one of the brightest OAVR star in 
the G C is a massive binary star (Mb ~ IOOMq). More recent 
works teafelski et al.ll200% found a few additional binary 0-star 
candidates, though no binaries had been found among the young 
B-stars. These studies either rely on detection of eclipsing binaries 
or through radial velocity detection of short-period binaries, and are 
therefore mostly sensitive to close binaries. It is therefore difficult 
to estimate the underlying binary fraction among the young stars 
in the Galactic center using these limited statistics. Taking them at 
face value would suggest a very high binary fraction, comparable 
or even larger than the high bi nary fraction obse rved in the filed for 
similarly young massive stars dSana et alj2012t) . 
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On the theoretical side, most binaries in the Galactic center 
are expected to be “soft” binaries and therefore have relatively 
short lifetime due to perturbations by other stars llBahcall & Wolf! 
I 977 I: Ozemoi & Dokuchaevlll982l : iHonmanll^OO^ Peretsll2009al : 


Alexander & Pfuh]||2014h .' or merge due to secular perturbations 


by the massive bla ck hole ( |Peretsll2009bl : lAntonini & Peretsll2012l : 
IProdan et alj|201^ . The binary fraction among the old stellar pop¬ 
ulation in the nuclear cluster is likely to be negligible, and we do 
not include it in any of our calculations. 

Binaries in the disk could, however, be “hard” in respect to 
the other disk stars, given the low-velocity dispersion in the disk. 
The young age of the disk also suggests that most of the initially 
formed binaries should still survive. Overall it is therefore likely 
that the binary fraction of the massive stars in the Galactic center 
stellar disk is high as in the field (> 50%), and that the binary 
fraction does not evolve much since its formation. 

As we discussed above, binaries are generally not very im¬ 
portant for the disk evolution, but we do consider them in some of 
our models for completeness, in which case we assume a binary 
fraction of 50%. 


3.3 Cusp and core models for the nuclear cluster 

A stellar cluster around a massive black hole is ex pected to evolve 
into a stellar cusp structure over a relaxation time jBahcall & Wolf! 
I1976L Il977h . For a spherically symmetric distribution of equal- 
mass stars analytic considerations and and N-body simulations 
show that in equilibrium the density profile of t he cluster has a 
powe r-law distribution, p oc r \ with 7 = 7/4 jBahcall & Wolll 
Il976h . Mass segregation in multi-mass clusters produce mass- 
dependent density profiles for the different mass stellar pop¬ 
ulations, with p{m) oc r where 7 (m) = 1.5 -f 

m/(4minax), with mm ax the mass of the heaviest stellar element 
jBahcall & Wolflll977h. 


Hopman & Alexandej j2006hh used Fokker-planck calcula¬ 


tions to show that 7 should be 1.4 for solar mass MS stars and 
up to 2 for stellar black holes (S BH). Similar result s from Monte- 
Carlo calculations were found bv iFreitag et alj j2006h However, ob¬ 
servations of red-giants in the Galactic center suggest a core-like 
structure in the inner regions of the GC (7 in the rang e 0 — 0.5; 
IDo et aljl2009l : iBartko et alj|201(il ; iGenzel et al.ll20ld Neverthe¬ 
less, it is not yet clear whether the distribution of red-giants re¬ 
flect the overall distribution of stars in the GC, and various mod- 
els had been suggested to explain both a “real” core distribution 
jMerritlll2010lt or an apparent one (i . e., only reflecting the distribu¬ 
tion of red-giants ; Dale et alll2009l ; lAmaro-Seoane & Cherj|201^ ; 
lAharon & Peretsll2015h. 

In light of the above discussion, we consider two possible 
models: (1) A cusp model dominated by a steep power-law dis¬ 
tribution of SBHs (7 = 2) and (2) A core model with a shallow 
power-law distribution (7 = 0.5). 


3.4 Stellar evolution and the long term evolution of disks 

In addition to our study of the role of the different NSC and disk 
components, we also considered the effects of stellar evolution, 
which become especially important once the long term evolution 
of the stellar disk is explored. We considered the evolution of disks 
with various properties (see below) both at times comparable to the 
observed stellar disk in the GC (up to 10 Myrs), as well as a longer 
term (100 Myrs) evolution. 


Table 1. The final mass and final type after Tms (of oach star) for different 
intervals of mass. However star with initial mass less than 6.5M0 has Tms 
greater than lOOMyr. 


Initial Mass [Mq] 

Final remnant mass (after Tms)IMq] 

Type of remnant 

30 <M ^ 120 

10 

SBH 

15 < M s: 30 

7 

SBH 

8 ^ M sg 15 

1.4 

NS 

6 sS M sj 8 

1 

WD 


Table 2. NSC Models. 


Model 

Star 

M{Mq) 

7 

PO (A70pc-^) 

CUSP I 

MS 

1 

2.0 

5.2 X lO'" 

CUSP II 

MS 

1 

1.4 

1.9 X 10® 


SBH 

10 

2.0 

1.8 X lO'^ 

CORE 

MS 

1 

0.5 

5.2 X lO'" 


We first considered simple models in which we did not intro¬ 
duce any stellar evolution, and assumed stars do not change over 
time. We then consider the effects of stellar evolution in a simpli¬ 
fied manner. For each stellar population in a given mass-bin we 
consider the appropriate MS lifetime for a star of such mass, Tms, 
according to stellar evolutionary models. After that period we re¬ 
place the star in the model with the stellar remnant it produces, 
assuming a simplified prescription as described in table [T] in par¬ 
ticular the continuous mass-loss process is simplified and the mass 
loss is assumed to be immediate. We only considered evolution of 
up to 100 Myrs and hence we do not consider the stellar evolution 
of stars with masses smaller than 6 M 0 . 


4 DETAILED MODELS 

In the following we describe the specific evolutionary models we 
discuss in detail. We consider a disk of stars orbiting a MBH with 
mass M, = 4 x 10®M 0 . The inner disk radius is = 0.05pc, 
and the external radius Rout = 0.15 pc (so AR — O.lpc). The 
initial velocity dispersion is assumed to be low (i.e. a cool thin disk 
with low eccentricity orbits) (initial velocity dispersion ajo <C Vk 
for any disk population j; we take than cryo = O.OSVk) . We take 
a disk s urface density of S (77) oc R~^ (following similar assump¬ 
tions as ICuadra et alJlrOOSl and IProdan et alj2015) . 

We consider seyeral models for the nuclear cluster, both cusp 
and core models. We study two different cusp model s. The CUSP 
I mod el is composed of single, solar mass stars (see IProdan et al.l 
I 2 OI 5 L for details) with a density profile profile 


p (r) = po 



1-f 



(17) 


with ro = 0.5 pc, po = 5.2 x 10® M 0 pc“® and 7 = 2 . 

The CUSP II model is a multi-component cusp, with two pop¬ 
ulations; MS stars of 1 M 0 , and SBHs of 10 M 0 . These have a 
power-law p oc r“^, distributions with 7 ms = 1-4 ysBH = 2 , re - 
spectively; following the results of lHopman & Alexandej j2006bh . 
Finally the CORE model is similar to CUSP I model but with 
7 = 0.5. The models properties, are summarized in Table|2 

We studied different models for various disk and cluster com¬ 
binations. A brief summary of the properties of the models is given 
in Table[3 including the disk mass function and its binary fraction. 
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Table 3. Summary of the properties of the disk and nuclear cluster models, including their mass functions and density profiles. 


# 

Disk 

Cusp/Core 

Masses in the cusp 

7 

Minin 

-Mmax 

r 

Binary fraction 

I 

O-stars with M* = 20 Mq 

- 

- 

- 

- 

- 

- 

- 

II 

O-stars with M* = 20 Mq 

CUSP I 

IMq ms stars 

1.5 

- 

- 

- 

- 

III 

O-stars with M* = 20Mq 

CORE 

IMq ms stars 

0.5 

- 

- 

- 

- 

IV 

Salpeter IMF 

- 

- 

- 

0.6Mq 

120Mq 

2.35 

- 

V 

Salpeter IMF 

CUSP I 

IMq ms stars 

1.5 

0.6Mq 

120Mq 

2.35 

- 

VI 

Salpeter IMF 

CORE 

IMq ms stars 

0.5 

0.6Mq 

120Mq 

2.35 

- 

VII 

Top-Heavy IMF 

- 

- 

- 

0.6Mq 

120Mq 

0.45 

- 

VIII 

Top-Heavy IMF 

CUSP I 

IMq ms stars 

1.5 

0.6Mq 

120Mq 

0.45 

- 

IX 

Top-Heavy IMF 

CORE 

IMq ms stars 

0.5 

0.6Mq 

120Mq 

0.45 

- 

X 

O-stars with M* = 20 Mq 

- 

- 

- 

- 

- 

- 

0.4 

XI 

O-stars with M* = 20 Mq 

CUSP II 

IMq ms /IOMq SBH 

1.4/2.0 

- 

- 

- 

- 


5 RESULTS 

In the following we present the evolution of stellar disks for the dif¬ 
ferent models considered. We first study simpler cases with single¬ 
mass disks and then progressively consider more and more realis¬ 
tic cases. These include disks composed of stars with a range of 
masses and different mass-functions, and the effects of different 
types of NSCs. In this case we also explore the differential evolu¬ 
tion of stars of different masses in the same disk and their stratifi¬ 
cation. We then consider the effects of stellar evolution on the disk, 
and the long term (100 Myrs) evolution of stellar disks. 

5.1 Evolution of Single-mass disks 

Fig-IUshows the evolution of the RMS eccentricity of the disk stars 
for disks with single-mass stars (models I, II, III and XI), for dif¬ 
ferent choices of mass (10 and 20 Mq stars). As expected the cusp 
heating has a stronger effect on lower mass stars. 

5.2 Dependence on stellar mass and the structure of the 
nuclear stellar cluster 

Fig. m shows the role of the cusp and its structure on the evolution 
of disk stars, and the role of the single-star masses. The final RMS 
eccentricity (after lOMyr) is shown for a range of single-mass disks 
(models I, II and III). Each point in the lines corresponds to a disk 
with the same total mass, but the mass of the stars composing the 
disk, indicated on the X-axis (same for all stars in the disk) differes. 
Disks composed of lighter stars can be seen to be more sensitive to 
the heating by NSC stars, with the cusp model being more effec¬ 
tive than the core one. In contrast, disks composed of heavy stars 
are cooled down by the NSC (compared with the model with no 
existing NSC). 


5.3 Multi-component NSCs 

In the previous section we studied the evolution of stellar disks and 
their interaction with single-mass NSC stellar population. In the 
following we consider a more complex and potentially more real¬ 
istic model for the NSC. In our multi-component NSC we follow 
model XI, where the NSC includes both SBHs and MS solar mass 
stars. Fig. [5] shows the disk evolution in such model. Though the 
number of SBHs in the cusp is relatively small, they can play an 
important role in heating a disk of massive 20 M© stars, where more 
massive SBHs lead to stronger heating of the disk. 

Fig. [6] shows that the inclusion of SBHs in addition to Solar 




Figure 3. Evolution of RMS eccentricity of single-mass disk stars. Upper 
panel shows models I, II, III and XL Bottom panel shows evolution of simi¬ 
lar models but with disk stars mass of IOMq . Also shown are the measured 
RMS eccentricity of the O-stars in the GC disk lYelda et alJ20l4l . 

mass stars in the cusp (model XI) leads to a more prominent disk 
heating, as expected. 

5.4 Evolution of multi-mass disks 

In models IV — IX we studied the evolution of multi-mass stellar 
populations in a disk, where we considered both disks with Salpeter 
IMFs as well as disks with top-heavy IMFs. For each model we 
divided the masses into 20 logarithmic mass bins for the Salpeter 
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Figure 7. Evolution of the RMS eccentricity of stars with different masses in the same disk. The upper panel corresponds to a disk with a Canonical IMF 
(Salpeter, F = 2.35) and the bottom one conesponds to a disk with a top-heavy IMF (F = 0.45). The different lines, ordered from bottom to top con'espond 
to the evolutionary times: 2, 4, 6, 8 and 10 Myrs, respectively. 




Figure 8. Evolution of the RMS eccentricity of all stars (of any mass) in each of the models. The left panel shows top-heavy IMF (F = 0.45) models and the 
right one shows models with a Salpeter IMF (F = 2.35). Solid lines correspond to models without any NSC contribution. 



Figure 4. The RMS eccentricity of disk stars after 10 Myr of evolution. 
The solid line corresponds to Model II, dot-dashed is for core Model III, 
and the dashed line corresponds to the case where no NSC exists. 


IMF cases; for better mass resolution we considered linear mass 
bins for the Top-Heavy IMF. For each mass bin we add a coupled 
equation ® , where we took the mass weighted average to represent 
the mass bin, with the appropriate number of stars in that mass bin. 

Fig. |7] shows that the disk evolution leads to a mass stratified 
structure of the disk, with lighter stellar populations being excited 
to higher eccentricities/inclinations (and larger scale height for their 
disk component) compared with populations of more massive stars, 
which are segregated to more circular orbits in the central part of 
the disk (lower scale height). Such stratification is also observed 
for models including the effect of the NSC, where only small dif¬ 
ferences exist as a function of the stellar masses; lower mass stars 
are more sensitive to the NSC heating (not shown), as expected. 

The top-heavy IMF models include a significantly larger num¬ 
ber of more massive stars, leading to a much faster heating of the 
stellar disk, as a l so ob served in the simple models considered by 
[Alexander et alj ( l2007h . As can be seen in Fig. [S] the NSC can 
further contribute to the disk heating, leading to ~ 10% higher 
RMS eccentricity, but its contribution becomes negligible in the 
top-heavy IMF cases, where the heating is dominated by the mas¬ 
sive stars in the disk. 
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Figure 5. Evolution of a stellar disk composed of 20 Mq stars, embedded 
in a multi-component cuspy NSC (model XI). The left panels shows the 
comparison between the disk-Hmulti-component cusp model and a model 
without an NSC component. The right panel shows the evolution of a disk 
embedded in cuspy NSCs, where different masses for the NSC SBH pop¬ 
ulation are cons idered. Also show n are the measured RMS eccentricity of 
the GC 0-stars iVelda et all20l4l . 

In order to compare these results with the measured RMS ec¬ 
centricity of the 0-stars in the GC we now consider only the RMS 
eccentricity of 3-40 Mq stars, more comparable to the 0/B stars 
observed in the disk today. For an age of 5 — 7 Myrs for the GC 
stellar disk we expect an RMS eccentricity of the order of 0.14, sig¬ 
nifica ntly lower than the measured (e) = 0.27 ± 0.07 jYelda et al.l 
|2014 . 


5.5 Stellar evolution and the long term evolution of stellar 
disks 

Fig.|T0]shows the long term (100 Myr) evolution of the stellar pop¬ 
ulations in disks with continuous (Salpeter or top heavy) IMFs, 
where simplified stellar evolutionary mass loss is considered (see 
section [J!4t . 

As each of the modeled disks evolve, progressively more mas¬ 
sive stars end their life and become lower mass stellar remnants. 
These (now) lower-mass remnants can then be heated much more 
effectively through mass-segregation processes. Once the more 
massive stars evolve the role of the cusp heating becomes more 
significant, as be seen from the comparison between the no-NSC 




Figure 6 . Evolution of RMS eccentricity for a disk of 0-Stars with multi- 
component cusp (model XI), each color is for different SBH mass. The 
dashed line is to compare to model I. On the left panel is plotted the evo¬ 
lution through the time. The right panel is the difference between the final 
eccentricity for each model compared to model I. 



Figure 9. Evolution of the RMS eccentricity of stars in the mass range 
3M0 < M < 40Mq corresponding to OB stars currently observed in the 
GC disk. 
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and cuspy NSC models. Moreover, with the absence of high mass 
stars intermediate mass populations begin to cool, when no cusp 
exists; but these same populations continue to heat up when a cusp 
population do exist. 


6 DISCUSSION 

In this work we studied the evolution of a stellar disk around 
a MBH. We followe d previous analytic models developed by 
[Alexander et al.l ( l2007h . which considered two-body relaxation and 
mass segregation processes, and extended the models to include the 
effect of binary heating, and role of the nuclear cluster interactions 
with the disk, using several plausible models for the NSC. We also 
studied realistic mass function models for the disk stellar popula¬ 
tions, and included simplified models for the stellar evolution of 
the stars. Finally, we also explored the long term (100 Myrs) evo¬ 
lution of such disks. In the following we discuss our results and 
their implications for the role played by the different processes and 
components in sculpting the disk evolution. 

6.1 Disk heating and mass stratification 

Similar to the results of [Alexander et al] ([2007h . we confirm that 
for isolated stellar disks, the disk heating is dominated by dynami¬ 
cal friction /mass-segregation processes in which the more massive 
stars heat the low mass stars. Correspondingly, the more massive 
stars are cooled by the interaction with the low-mass stars, but two- 
body relaxation of the massive stars by themselves still keeps heat¬ 
ing them. Models with top-heavy mass function include a larger 
fraction of massive stars, and therefore allow for more rapid heat¬ 
ing of the disk, producing larger RMS eccentricities. 

Our models considered detailed mass-functions for the stars 
in the disk, allowing us to estimate the evolution of different stellar 
populations. As seen in Fig. [7] the disks develop mass stratification 
where higher mass stellar populations are expected to have signif¬ 
icantly lower disk height scale and more circular orbits compared 
with the low mass stars in the disk. The existence, or lack of such 
mass stratification in the GC stellar disk could therefore be an im¬ 
portant tool in assessing the processes involved in the evolution of 
the GC disk and its initial conditions. 

6.2 The role of the cusp 

The effect of the NSC on the evolution of the stellar disk was not 
considers before in analytic models. The NSC stellar population 
has a high velocity dispersion, i.e. it is an effectively “hotter” pop¬ 
ulation than the “cool” stellar disk population. The NSC can there¬ 
fore potentially heat the stellar disk, even when it is composed of 
low-mass stars. 

Our results suggest that a dense cuspy NSC can play a non- 
negligible, though limited role in heating the stellar disk. When the 
disk contains a large population of massive stars, these dominate 
the disk heating and evolution, and the NSC plays a relatively mi¬ 
nor role. Nevertheless, even young disks in which massive stars still 
exist are affected by the NSC, with top-heavy IMF disks being rela¬ 
tively little affects, but Salpeter-IMF disks showing non-negligible 
heating up to ~ 10% higher RMS eccentricities. 

The NSC effects are more prominent in cases where the disk 
is dominated by low mass stars. These include unrealistic cases of 
disks composed of low single-mass stars, but they also relevant for 
the long term evolution of stellar disks with realistic IMFs. In the 


latter case stellar evolution leads to the transformation of massive 
stars into low mass stellar remnants over time. The effective mass- 
function of such evolving stellar disks therefore becomes progres¬ 
sively centered around lower mass stars and stellar remnants which 
replace the previously existing massive stars. Correspondingly, the 
NSC heating of the cusp plays a progressively more dominant role 
in the disk evolution as can be seen in Fig.llOl 


6.3 The role of binary heating 

In this study we extended the analytic study of stellar disk evolu¬ 
tion to include binary heating processes. In principle the binding 
energy stored in binaries is significant; however, the rate of energy 
exchange between the binaries and the single stars in the disk is 
slow. Therefore, binary heating is inefficient in heating the stel¬ 
lar disk, compared with two-body relaxation and mass-segregation 
processes. These results are consisten t with N-body simul ation re¬ 
sults which included binaries done bv [Cuadra et al.l ( [2008h . 

Note that we consider the overall averaged evolution of 
the stellar disk population; atypical single-binary encounters with 
short-period massive binaries can lead to strong kicks ejecting stars 
at high velocities, ev en beyond the escape velocity from the NSC 
jPerets & Subi[[20l^ . Such encounters may produce a small num¬ 
ber of outlying stars at highly eccentric/inclined orbits, but are not 
likely to significantly affect t he overall evolutio n of the disk and its 
averaged properties (see also [Perets et alj2008ll . 


6.4 The long term evolution of steallar disks and their 
signature 

We find that the disk evolution over times as long as 100 Myrs is 
slow. Two-body relaxation is too inefficient for the disk to assimi¬ 
late into the nuclear cluster on such timescales. The disk structure 
is expected to keep its coherency, and be observed as a relatively 
thin disk even at 100 Myrs. This suggests that the existence of 
older disks formed before the one currently observed in the GC 
might still be inferred from the stellar kinematics of older, lower 
mass stars, which may still show a disk-like struture, unless de¬ 
stroyed/smeared by other non-two-body relaxation processes. 


6.5 Other physical processes (and caveats) 

Our analytic approach considered the role of two-body relaxation 
and binary heating, and included realistic disk and NSC compo¬ 
nents. However, secular processes and possible collective effects 
are beyond the scope of such a model; we briefly review studies of 
these processes. _ 

Resonant relaxation processes ([Rauch & Tremain3 [ 19961) 
could potentially lead to fast evolution of disk stars inclina¬ 
tions through the vector resonant relaxation (and slower evolu- 
tion of the eccentricities), espe ci ally in the inner par t s of th e disk 
( [Hopman & AlexandeJ [2006a[ ). [Kocsis & Tremain3 ( [201l[) sug¬ 
gested that resonant relaxation can produce significant effects, how¬ 
ever they assumed an unrealistic effective mass of NSC stars, which 
is 5-10 times larger than expected in the GC (even when including 
stellar black holes). Further study of this effect under realistic con¬ 
ditions could shed more light on this issue. 

Fast two-body relaxation as well as secular evolution of 
disk stars could also be i nduced by massiv e coherent compo¬ 
nents, massive perturbers ( [Perets et alj [TOOtI) such as molecular 
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Figure 10. The long term evolution of the RMS eccentricity of disk stars in over 100 Myrs. Simplified stellar evolution is included (see main text). Top panels 
show the evolution of a disk with a a Top-Heavy IMF (F = 0.45). Bottom panels show the evolution of a disk with a Salpeter IMF (F = 0.45). The left and 
right panels correspond to models without an NSC and with a cuspy NSC, respectively. The lines from bottom to top (initially) correspond to increasingly less 
massive stars. The solid, dashed, dashed-dotted and dotted lines correspond to star that evolve to become SBHs, NSs, WDs or MSs, respectively. 


clouds, s tellar clusters, intermediate black hole or a second stel¬ 
lar disk iLdckm ann et alJboogI : I Yu et al.ll2007l : iHaas et ^1201 ll : 
iMapelli et al.l2013l) . Whether such objects exist/had existed in such 
configurations as to influence the evolution of the stellar disk is still 
actively studied. 

iMadigan et al.l ( l2009l) considered the collective eccentric disk 
instability effect. Such effect could lead to rapid change in the dis¬ 
tribution of eccentricit ies in the disk, but it req uires initially highly 
eccentric disk (see also lGualandris et alj2012h . Such a process may 
give rise to a significant population of highly e ccentric stars in th e 
disk; which are not seen in most recent studies jYelda et al.ll20l4 . 

It is important to note that by their nature, most of the collec¬ 
tive and secular effects mentioned above affect both low mass and 
high mass stars in a similar manner, leading to similar eccentricity 
distribution. Discerning the kinematic properties of low mass vs. 
high mass stars is therefore highly important in order to understand 
the roles of such processes compared with stellar two-body relax¬ 
ation processes that lead to mass stratification, as shown here in 
details. In particular, non-dependence of the stellar kinematics on 
the stellar masses would strongly suggest collective effects domi¬ 
nate the stellar disk evolution. 


6.6 Comparison with observations of the Galactic center 


The observed eccen tricity of 0-stars in the stellar disk in the GC, 
(e) = 0.27 ± 0.07 dYelda et al.ll2014l) . is higher than that obtained 
in the realistic models considered here. This suggest either initial 
conditions with significantly higher eccentricities for the disk stars, 
or that other processes besides two-body relaxation play an impor¬ 
tant role in the disk evolution. 

Detailed studies of the eccentricity/inclination distribution of 
disk stars as a functi on of mass is not currently available, (but see 
iMadigan et al.ll2Ql^ for important global statistical trends). More¬ 
over, lower mass B-stars in these regions may have a different ori¬ 
gin than th e disk, and may ha ve been captured following a binary 
disruptionj|Perets_et^ 200^j_siKh stars would have high eccen¬ 
tricities dPerets & Gualandrisll2010h and may mask the real distri¬ 
bution of the disk stars eccentricities. Nevertheless, a focused study 
on stars most likely related to the disk, and their eccentricity-mass 
relation is highly desirable to asses the existence of possible mass 
stratification in the stellar disk. 

Finally, two body relaxation processes suggest a factor of two 
ratio between the eccentricities and inclination of stars in the disk, 
while different relation might be expected from other processes 
(e.g. secular Kozai-like evolution induced by an inclined massive 
structure). A good handle on the eccentricity vs. inclination from 
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observation could give additional clues on the dynamical processes 
involved. 


7 SUMMARY 

In this study we analyzed the evolution of a stellar disk around 
a massive black hole due to two-body relaxation processes and 
binary heating. We explored realistic mass functions for the disk 
stars, and included the effects of two-body relaxation by the nu¬ 
clear cluster stellar population in which the disk is embedded. We 
also considered the effects of stellar evolution and studied the long 
term evolution of such disks. The disk evolution is dominated by 
dynamical friction from the massive stars in the disk, while binary 
heating plays only a negligible role in the disk evolution. The nu¬ 
clear cluster plays a minor role in the disk evolution as long as a 
large population of massive stars exist, as in the case of a young 
disk with a top-heavy mass function; it can play a more significant, 
though still modest role in heating a disk with a Salpeter mass func¬ 
tion. Due to stellar evolution the effective mass-function of the disk 
becomes progressively centered around lower mass stars and stel¬ 
lar remnants which replace the previously existing massive stars. 
Correspondingly, disk heating by the nuclear cluster plays a pro¬ 
gressively more dominant role in the disk evolution. 

We find thaf significant mass stratification arises from mass- 
segregation processes in the disk. This could serve as a signature 
for two-body relaxation dominating the disk evolution. In contrast, 
collective and secular effects are typically insensitive to the masses 
of individual stars in the disk, and would not produce mass stratifi¬ 
cation. 

The observed RMS eccentricity of the O-stars in the stellar 
disk in the Galactic center (~ 0.27) is larger than obtained in any of 
our models with realistic conditions, suggesting either an initially 
hot and/or eccentric disk or that other secular/collective effects play 
an important role in the disk evolution. 
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